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Abstract 

We introduce conformal multi-matrix models (CMM) as an alternative to conven- 
tional multi-matrix model description of two-dimensional gravity interacting with 
c < 1 matter. We define CMM as solutions to (discrete) extended Virasoro con- 
straints. We argue that the so defined alternatives of multi-matrix models represent 
the same universality classes in continuum Umit, while at the discrete level they 
provide explicit solutions to the multi-component KP hierarchy and by definition 
satisfy the discrete VF-constraints. We prove that discrete CMM coincide with the 
{p, g)-series of 2d gravity models in a well-defined continuum limit, thus demonstrat- 
ing that they provide a proper generalization of Hermitian one-matrix model. 
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1 Introduction 



Matrix models are used nowadays to describe non-perturbative partition functions, which 
interpolate between various sets of two-dimensional conformal models coupled to 2d grav- 
ity, and thus serve to approach the old-standing goal of devising a universal partition 
function of entire string theory. Considerable progress so far is reached in unification 
of all tachyon-free string models, associated with c < 1 minimal conformal models, i.e. 
the models from {p, g)-series plus 2d gravity. Namely, in it has been suggested that 
the non-perturbative partition function Zp[Tk], interpolating between all the models with 
given pJ3 possesses the following properties: 

(i) Zp[Tk] is a r-function of p-reduced KP hierarchy, which is completely independent 
of all Tfc with A; = mod p; 

(ii) Zp[Tk] satisfies an infinite set of differential equations of the form 

Wt^Zp = 0, a = 2,...,p; i>l-a, (1) 

where wf'^'* are (harmonics of) the generators of Zamolodchikov's Wp-algebra ^, (see 
sect. 4 below for explicit expressions of wf"^ in terms of T and d/dT). These two state- 
ments (and actually the second one alone) are enough to define Zp[T] completely. Further 
interpolation between different values of p is described in terms of Generalized Kontsevich 
Model (P, 1^ and references therein). 

Though the above formulated suggestion of can be just used as a reasonable defini- 
tion of non-perturbative partition functions (and these are the corollaries of this definition 
which should be studied to proceed further with the string program), there are still many 
interesting but unresolved problems, which concern the original motivation of that sug- 
gestion. The most direct motivation is usually supposed to come from the study of the 
multi-scaling continuum limit of discrete multi-matrix models (the number of matrices be- 
ing p — 1). While being completely understood in the simplest case of p = 2 (one- matrix 

model) 1]^, such continuum limit was never honestly studied for any p > 2. This prob- 
^ Partition function of a given (p, q) string model (i.e. the generating functional of all the corre- 
lation functions in the model) arises from Zp[Tk], if all the time- variables — 0, except for Ti and 
Tp^q. Zp,q\T] ^ Zp[Ti, 0, . . . , 0, Tp^q, 0, . . .]. Existing formalism does not respect explicitly the symmetry 
between p and q. 
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lem attracted certain attention, but remains unresolved because of considerable technical 
difficulties. 

In the present paper we are going to argue that there is a way around all these 
difficulties. Moreover, we continue to advocate the central idea of 0] that one has to 
look for the origin of all essential features of continuum double-scaling limit of (multi-) 
matrix models in their discrete counterparts. From this point of view one accepts that the 
conventional multi-matrix models, as defined in |P, ^ |ll| are not the simplest possible 
discrete representatives of the relevant classes of universality. This is more or less clear 
already from the fact that these do not satisfy any simple Ward identities at the discrete 

- is 



level for at least, the structure of these identities — so called VT-constraints 12 



rather different from the ly-constraints (|1])). Moreover, the advantage of being solved by 
orthogonal polynomial technique is not enough to present them in the form which respects 
the whole integrable structure. 

Below we shall con&m the suggestion of |l| that much better representatives of the 
same universality classes are provided by somewhat different "multi-matrix" models or 
just solutions to integrable hierarchies and string equations which can be written in the 
multiple-integral form essentially different from conventional multiple integrals of multi- 
matrix models (by "interaction" of the nearest neighbours as well as by specific choice of 
integration contours). We shall call them conformal multi-matrix models (CMM), because 
of their relation to the formalism of conformal field theory. These models are by definition 
constructed to satisfy H^-constraints at discrete level. Moreover, some especially inter- 
esting solution to these discrete VT-constraints immediately at the discrete level possess 
a rich integrable structure which is in some sense less rough than the integrable structure 
found in the case of ordinary multi-matrix models [§, More concretely, the partition 
function of the standard multi-matrix models is a r-function of Toda lattice hierarchy 
in the ffist and the last times, the other times describing only the parameterization of 
the point of the Grassmannian. Thus, this case corresponds to extremely non-economic 
usage of the whole variety of parameters (times) in the problem. The way to rule out 
this drawback is to involve all these parameters as new times of a more rich hierarchy, 
namely, some multi-component hierarchy. It turns out that our CMM are related to the 
multi-component KP hierarchy of the SL{p) type (or generalized AKNS type 0)- 
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Thus, we can claim that the proper viewpoint on the one-matrix model is to consider it 
reduced 2-component KP hierarchy. 
In what follows we shall discuss some generic (especially integrability) properties of 
CMM in more details than it was done in and give a detailed description of the 
relevant continuum limit (analogous to the presentation of [0] for p = 2, i.e. including the 
definition of Kazakov variables tk ^ Tk, reduction, rescaling of partition function and the 
transformation of discrete VT-constraints into the continuum ones (|l]), see also |ll6[)- As 
a result, we argue that the continuum limit 

i/p 

<i zz'Ar''""'\ti.\ ] 

is exactly a solution to eqs.(p, i.e. the relevant partition function of the c < 1 string 
theory. 

The plan of the paper is as follows. In the sect. 2 we remind (and make more explicit) 
the construction of CMM from ref. and discuss the form of (Zamolodchikov's) W- 
constraints in these discrete models. In the sect. 3 integrability properties of 
are discussed. The sect. 4 is describes the p- reduction of CMM, the analogue of Kazakov 
variables and the transformation of discrete into continuum ly-constraints mainly along 



hm {Z™-^N}'/^ = Z,[T,] (2) 



lines of the ref.[|TB|. Concluding remarks are given in the sect. 5. 



2 Conformal mult i- matrix models 

In this section we shall remind the ideas of and give the definition of CMM. First, 
we show that the simplest example of discrete Hermitian 1-matrix model can be easily 
reformulated in these terms. 

Indeed, Hermitian one-matrix model (p = 2) can be defined as a solution to discrete 
Virasoro constraints: 

-^n^2,AfM = 0, n > -1 

oo 

L„ = ^ -I- dVdtadh (3) 

k=0 a+b=n 
dZ2,N/dto = -NZ2,N 

The Virasoro generators (|^) have the well-known form of the Virasoro operators in the 
theory of one free scalar field. If we look for such solution in terms of holomorphic 
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components of the scalar field 



[Jn,Jm] =n5n+m,0, [^,P] = 1 (4) 

the procedure is as follows. Define vacuum states 

Jk\0) =0, {N\J_k = 0, k>0 

p\0) =0, {N\p = N{N\ (5) 

the stress-tensor0 

T{z) = lmz)]' = J2T^z-^-\ T„ = 1 ^ J_, J,^„ + 1 ^ J, J,, (6) 

^ ^ fc>0 '^ + b=n 

a,b>0 

Tn\0) = 0, n>-l (7) 

and the Hamiltonian 

Hit) =^Y.tkJk=i V{z)j{z) 

V{z) =Y.tkz\ 3{z) = ^dct^{z). (8) 

fc>0 

Now one can easily construct a "conformal field theory" solution to (j^) in two steps. First, 

L„(iV|e^W... = (iV|e^WT„... (9) 

can be checked explicitly. As an immediate consequence, any correlator of the form 

(iVle^^'^GlO) (10) 

(iV counts the number of zero modes of G) gives a solution to (|^) provided 

[T„,G]=0, n>-l (11) 

Second, the conformal solution to ( [TT| ) (and therefore to (|^)) comes from the properties 
of 2d conformal algebra. Indeed, any solution to 

[T(z),G] = (12) 



^For the sake of brevity, wc omit the sign of normal ordering in the evident places, say, in the expression 
for T and W in terms of free fields. 



is a solution to ([TT|), and it is well-known that the solution to (|T^) is a function of screening 
charges 

Q± = fj± = fe^^t (13) 

With a selection rule on zero mode it gives 

G = exp Q+ ^ l^Q^ (14) 

(Of course, the general case must be G ~ Q^^^Q'^ but the special prescription for 
integration contours, proposed in [|l|, implies that the dependence of M can be irrelevant 
and one can just put M = 0.) In this case the solution 

Z2,^[t] = (iV|e^WexpQ+|0) (15) 

after computation of the free theory correlator, analytic continuation of the integration 
contour (see more detailed discussion of this point below) gives well-known result 

r ^ 

Z2,N = l[dziexp{-Y.hz'[)A%iz) = 

1=1 

= {NNo\ U{N))-^ J DM exp (- ^ t^M^^) (16) 

TV 

Aiv = Ylizi - Zj) 

in the form of multiple integral over spectral parameter or integration over Hermitian 
matrices. 

In the case of p = 2 (Virasoro) constraints this is just a useful reformulation of the 
Hermitian 1-matrix model. However, in what follows we are going to use this point of 
view as a constructive one. Indeed, instead of considering a special direct multi-matrix 
generalization of ( P^ ) |rT| one can use powerful tools of conformal theories, where it 



is well known how to generalize almost all the steps of above construction: first, instead 
of looking for a solution to Virasoro constraints one can impose extended Virasoro or W- 
constraints on the partition function. In such case one would get Hamiltonians in terms 
of mu/iz-scalar field theory, and the second step is generalized directly using screening 
charges for W^-algebras. The general scheme looks as follows 

(i) Consider Hamiltonian as a linear combination of the Cartan currents of a level one 
Kac-Moody algebra Q 

i7(t«,...,t(-°'^^))= t^^^.J,, (17) 

A,fc>0 



where {n^} are basis vectors in Cartan hyperplane, which, say for SL{p) case are chosen 
to satisfy 

1 ^ 
P j=i 

(ii) The action of differential operators Wf"'* with respect to times {t^^^} can be now 
defined from the relation 

|^W(^|e^^(W)... = (iV|e^«*»w;'^)... , a = 2,...,p; z > 1 - a, (18) 

where 

W(^\z) =^[/XA9c^(z)r + ... (19) 

A 

are spin-a W-generators of Wp-algebra written in terms of rank ^-component scalar fields 
§• 

(iii) The conformal solution to (||) arises in the form 

= (ATle^^W^GlQ^^^llO) (20) 
where G is an exponential function of screenings of level one Kac- Moody algebra 

Qio^) = j ji.-) = j e'^4> (21) 

{ck} being roots of finite-dimensional simply laced Lie algebra Q. (For the case of non- 



simply laced case see [|T3- Below Q = SL{p) if not stated otherwise.) The correlator (^ 



is still a free-field correlator and the computation gives it again in a multiple integral form 

a [i=l y A,fc>0 

X n nn(4"^-^f)"^ (22) 

(a,/3) i=lj=l 

The expression ([2^ ) is what we shall study in this paper: namely the solution to discrete 
VT-constraints (|l]) which can be written as multiple integral over spectral parameters 
{-^i"^} (this integral is sometimes called "eigenvalue model"). The difference with the 
one-matrix case (|I6|) is that the expressions (^2]) have rather complicated representation 



in terms of multi-matrix integrals. Namely, the only non-trivial (Van-der-Monde) factor 
can be rewritten in the (invariant) matrix form: 



i=ij=i 



det{M(") O/-/® M^'^)} 



(23) 



where / is the unit matrix. Still this is a model with a chain of matrices and with closest 
neibour interactions only (in the case of SL{p)). 

The purpose of this paper is to show that CMM, defined by (^2]) as a solution to 
the VT-constraints has indeed a very rich integrable structure and possesses a natural 
continuum limit. To pay for these advantages one should accept a slightly less elegant 
matrix integral with the entries like (p^). 

The first non-trivial example (which we use as a demonstration making all our state- 
ments in the paper more clear) is the p = 3 associated with Zamolodchikov's PVs-algebra 
and serves as alternative to 2-matrix model. In this particular case one obtains 

1 

' fc>0 



V 2 uyn 



(24) 



Wj,^^ = Ln = Y.{ktkd/dtk+n + khd/dh+n) + 



fc=0 



a+b=n 



(25) 



Wi^^ = E {ktkltid I dtk+n+i - ktkltid/dtk+n+i - 2ktkltid / dtk+n+i) + 

k,l>0 



+2E 



A;>0 



E {ktkdydtadU - ktkd^/diadh - 2khd^/dtadh) 

a+b=n+k 
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+ 3 E {d^'/dtadm.-d^/dtadm,), (26) 

a+b+c=n 

where times tk and correspond to the two orthogonal directions in S'L(3) Cartan plane. 
(We use the standard specification of the Cartan basis: e = q;i/-\/2, e = -\/3i^2 /"\/2- This 
basis is convenient for discussing integrability properties of the model in the sect. 3, for 
the continuum limit we will use another basis in Cartan plane connected with t ± t.) In 
this case one has six screening charges Q'-^"*-' (i = 1, 2, 3) which commute with 



W(2)(z)=T(z) = -[90(;.)P 



(27) 
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and 

W(^)(^) = ^(MA50(z))^ (28) 

A=l 

where fi^ vectors of one of the fundamental representations (3 or 3) of SL{3). 

The particular form of integral representation (^) depends on particular screening 
insertions to the correlator (pOl). Following ref. we will concentrate on the solutions 
which have no denominators. One of the reasons of such choice is that these solutions 
possess the most simple integrable structure of the forced hierarchy type |1^ (see 
sect. 3), though the other ones can still be analyzed in the same manner. 

One of the possible ways to avoid the problem with other solutions is the following 

prescription with the integration contours. First, as in any conformal theory, one has 

to consider closed contours in the definition of screening charges (|T3p, (^TJ). Then, the 

denominators can be easily integrated out using the ordinary residue technique. However, 

in the final result, it is necessary to continue contours to their specific locations defining 

any concrete model. Say, in the simplest example of the Hermitian one-matrix model the 

integration runs along real line. This is not too simple procedure, because only after this 

the solution starts to be non-trivial. (In a sense, with properly defined contour Q^^^\0) = 

/ jW|0) ^ |0), what would be the case for closed contour). Thus, the prescription to 
J R 

get rid of denominators implies that one first integrates over some part of the spectral 
variables and then uses an analytic continuation. However, even besides this prescription 
the integrable structure of generic solutions (^) can be analyzed in the free-fermion 
formalism and we shall return to this problem in the sect. 3. 

The simplest solutions which have no denominators correspond to specific correlators 

^Sr[W] = (iV|e^«*»nexpg.jO) (29) 

i 

when we take to be "neighbour" (not simple!) roots: (ckjCKj) = 1. In the case of SL{3) 
this corresponds, say, to insertions of only Qai and 

Z^^'},[t,t\ ^ Zm,nM = ^(iV,M|e^(*'*")(Q(-))^(g-)*^|0) = 

1^ J l[dxdy^exp {-Y.[Vix.) + Viy.^]) A^(x)Ai,(y) Hix, - y,) (30) 



N\M\ 

III dT.:dlh P">fn I 

N\M\ 

This expression can be examined for its integrable properties by methods developed in 
]T8[ (see sect.3). 



Before we proceed to a detailed investigation of the integrable structure of CMM, let us 
make two remarks about possible generalizations. The first one concerns supersymmetric 



matrix models. This problem was first examined in where the authors looked for a 
solution to the system of equations L^Z = and GmZ = 0, with generators {L„, Gm} 
forming the = 1 superconformal algebra. The solution was found in the form of a 
multiple integral over even and odd (Grassmannian) parameters, in the form similar to 
that considered above. It is necessary to point out that in our language this is nothing but 
trivial generalization of the one-field case which has to be substituted by scalar superfield. 
Then the insertion of screenings of = 1 superconformal algebra immediately leads 
to the result by ||20[. From this point of view the real problem with supersymmetric 
generalization can arise only in the N = 2 case because of the lack of appropriate screening 
operators. 

The second possible generalization is to the case of "deformed" Virasoro— con- 
straints via Feigin-Fuchs-Dotsenko-Fateev procedure. In such case formulas (^) and (pO]) 
are still valid, one has only change the definitions (|]), (|19|) to 

T = ^[d(l){z)]' + aod'<P (31) 

(and corresponding formulas for higher W-algebras with c < rank^) with correspond- 
ing deformations of generators (|I]) written in terms of differential operators. Deformed 
generators will now commute with new screening charges 

Q± = f e"±*, a± = ao± (ao + 2)^/' (32) 

in the case of deformed p = 2 theory (minimal models). The only problem which might 
arise here is connected with the integrability of such system — expressions (|3^) have no 
natural representation via free fermions for generic a^. To avoid this problem one should 



consider a j-differential system |2T|. After the Miwa transform |^2[ of the time- variables 



^'« = rE«.e7' (33) 



our correlators acquire the form of 
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with 

^i{z) = exp[ai(j){z)], Ai = a-/2 - 0^00 (35) 

One can choose {aj} in such a way that {^j{z)} become primaries of minimal models, in 
such case the correlators will satisfy certain differential equations following from the 
null-vector conditions (in the particular case of 2-level degeneration they will be identical 
with the Virasoro constraints in Miwa variables). This can be also easily generalized to the 
W-case of several scalar fields and several sets of times. In this case one should introduce 
in (^) vectors and aj. This leads to independent differential equations which again 
correspond to the independent null-vector conditions. Such construction trivially reveals 
the origin of the correspondence between Virasoro constrained r-functions and correlators 
in the minimal models. 

Unfortunately the construction is valid in its present form only for discrete models 
(solution to discrete constraints) and is rather difficult to generalize to continuum models. 
The reason is that in the latter case we have to impose the twisted boundary conditions 
on scalar fields and thus (i) there is nothing like selection rule (0) which in the discrete 
case is provided by a zero mode and (ii) eq.(|^) is no longer satisfied. We are going to 
return to this problem in a separate publication (see also p^). 



3 Integrability of conformal mult i- matrix models 

In this section we are going to present a detailed investigation of integrable structure of 



CMM, based on the application of the general technique of ||T8[. In the subsect.3.1 we will 
derive a determinant representation for CMM partition function which is an indication 
to the fact that we deal with an integrable system and the partition function is a certain 
r-function. As a direct result of the determinant formula we obtain the Hirota bilinear 
relation satisfied by partition function of CMM, which is a generalization of that one for 
the Hermitian one-matrix model. In the subsect.3.2 we are going to apply the fermionic 
operator formalism for CMM which brings us to a conclusion that the partition function 
of CMM is a r-function of a multi-component Kadomtsev-Petviashvili hierarchy which 
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obeys the constraint 

E9/ae)r(^)({t}) = 0, n = l,2,... , (36) 

k=l 

where p is the number of components. This r-function satisfies l^-constraints (|l|) and is 
a discrete counterpart of the analogous statement in the continuum hmit found in . 

3.1 Determinant representation 

Let us first remind briefiy the resuhs of for the simplest example of p = 2. The 
partition function of p = 2 case ([l6| ) can be written in the form 



Z,Mt) = detJd^^^-'C{t)] = Mt) (37) 

with 

d,„C{t)=dl[C{t) (38) 

which implies that it is a r-function of the Toda chain hierarchy (or a r-function of AKNS- 
reduction of 2-component KP hierarchy see ref. and sect. 3. 2 below for details). 



Eq.(^) means that C{t) just has an integral representation 

C(t) = J dfi{z)expY.tkz'', (39) 

where is some measure; these are the Virasoro constraints which fix the concrete 
measure and the contour of integration in (|39D . The determinant form (|37D is an ex- 
plicit manifestation of the fact that the partition function does satisfy the Hirota bilinear 
relations, the simplest one of which in this particular case takes the form 

= rN+iit)r^-iit) (40) 

The fact that in eq.(|37|) we are dealing with a (forced) Toda chain reduction of generic 
Toda lattice (two-component KP hierarchy) is refiected in the specific feature of the matrix 

in (13): Q,=a+j MMM- 



Now we are going to generalize ( PTj ) and (^Op. In the case of the p = 3 model ( POD one 
has to introduce two functions instead of (Bol): 



C{t) = Jdz exp[-V{z)], C{t) = Jdz exp[-V{z)] (41) 
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where 



and 



k>0 k>0 



The determinant representation has now the form {d = d/dti, d = d/dii] 



(42) 



det 



C 

dC 



dC 



d^-^C C 



dc 



dc 



qN+M-Iq qN+M(j q2N+M-2q qN+M-1^ qN+Mq q2N+M-2(J 

= TN,Mit,t) (43) 



which is exactly the double- Wronskian representation of a r-function p6[] . 

We can prove the determinant representation (^) as follows. Instead of eq.(|5D|) let 
us consider shghtly more general expression by introducing the external sources {Pi} and 

m-. 

1 



Z^T^[t,t\f3,(3] 



N\M\ 



X 



j Wdxidyie^^ {'YJ^^^i) ' + ^ Ajv(x)AAf(y)A(x, y), (44) 



where A(x,?/) = /\N{x)^M{y) Y{i,j{xi-yj) is (A^ + M) x {N + M) Van-der-Monde deter- 
minant. Using the derivatives with respect to (3i and Pi one can get rid of A{x,y) in the 
integrand thus obtaining 

^3™ A^] = Aidp,d^)Fit,P)F{iJ), (45) 

where F and F have the specific form of Kontsevich-like integrals : 

Fit,P) = J nc^a;,exp [-Y,[V{xi) - PiXi]) A^ix) , 

F{IP) = J n^a;,exp (^-Y^[V{xi)-P,x,]) A^{x) . 

Using the trick with the differentiation over /? and f3 again one can represent these ex- 
pressions in the determinant form (compare with ^): 

F{t,P) = det[di-'C{t,Pj)] = det[dl;'C{t,P,)], z,j = l,...,N, (46) 
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F{t,P) = det[d^^^'C{tJj)] = det[dl;'C{t,Pj)], t,j = l,...,M, (47) 

where 

C{t,(3) = j dz e^^[-V{z) + I3zl C{tJ) = J dz exp[-V{z) + (3z]. 

Substitution of eqs. ( ^U]) and ( |^ ) into eq. (^31) leads to the expression which contains the 
sum of MINI terms analogous to (|43|) (each term corresponds to a particular permutation 
of {Pi} and {Pi}). Finally, in the limit of Pi = Pi = it reproduces the eq. (^31) . 

Now the generalization of (^) for p > 3 is quite obvious: r-function has "(p — 1)- 
tuple" Wronskian form for {J2 Ni) x [J2 Ni) matrix (multiplied by the factor H^J) with 
corresponding = / exp|-l'.(.)ld. (, = 1, . . . ,p - 1). 

From representation (^Sf ) it is easy to derive that 



J- log Tn^m [t, t) = 2 — TTTa 



dtidti ' ' 'T'NM^t.i) 

This expression can be also easily extended to the (p — l)-matrix case: 

, /, T{..,N,+l,-,Nj~l,..}{t,t)r{..,N,-l,..,Nj+l,..}{t,t) 

■log nN^}{'t,t) = 2 — 77T\ 



This equations of motion can be immediately derived from the corresponding determi- 
nant representation like the case of Wronskian solutions of KP- or Toda lattice hierarchies 
(see, for example, [^]) and are the first Hirota bilinear equations generalizing (^Op. 

3.2 Fermionic representation 

Now we shall proceed to the representation of the solutions to CMM in terms of free 
fermion correlation functions. Such a representation (invented for integrable hierarchies 
by Kyoto school allows one to establish some of the properties of the system under 
consideration in an elegant and "physical" way. 

Again, first we are going to show that the solution to Hermitian one-matrix model is 



nothing but AKNS-reduction of 2-component KP hierarchy |T8| . Indeed, the r-function 
of 2-component KP hierarchy is by definition the correlator 

ri^U^, y) = {N, Mle^'^-^y^GlN + M, 0) (50) 
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where 

/f(x,y) = 5:(x,j«+i/,jf ) (51) 



fc>0 



= J2 Jk^z-^'^ = (52) 

ij^'\z)ij^'>{z') = ^^ + ... . (53) 

Now we are going to demonstrate that (pTSf ) is equivalent to (^) for certain G for which 
(50) depends only on the differences Xk — Vk- To do this we have to make use of the 
free-fermion representation of SL{2)k=i Kac-Moody algebra: 

Jo = i(v;(^V^'^* - v^^'V^'^*) = ^(^^'^ - -/^'^) 

J+ = ^(2)^(1)* J_=^(iV(')* (54) 



Now let us take G to be the following exponent of a quadratic form 

G = : exp (^J ^(^V^'^*) : (55) 



The only term which contributes into the correlator (pD| ) due to the charge conservation 
rule is: 

Gn,m = Gn,-nSm,-n = ]yy • (y V-^'V^'^*] : hl,-N (56) 
Now we bosonize the fermions 



and compute the correlator 

TS\x,y) ^ T^tAx,y) = ^(iV,-iV|exp (xixkA'^ + ykA'^)] (/ \0) 



= ^(iV,-iV|exp(/[X(z)J«(z) + r(z)j(2)(^)]) (/ : exp(0i - 02) :) |0) 
Introducing the linear combinations -\/20 = 0i — 02, "\/20 = 0i + 02 we finally get 



T^\x,y) = ^(exp (-^f[X{z) + Yiz)mz)\) x 



X 



(iV|exp('i=^[X(z)-r(^)]a0(z)') (I ■.expV2<p-y \0) = Tl^\x-y) (5^ 
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since the first correlator is in fact independent of x and y. Thus, we proved that the 
r-function ( [50| ) indeed depends only on the difference of two sets of times {x^ — Vk}- So, 
we obtained here a particular case of the 2-component KP hierarchy (^) and 

(i) requiring the elements of Grassmannian to be of the form (|55|) we actually performed 
a reduction to the 1-component caseQ; 

(ii) we proved in (|58|) that this is an AKNS-type reduction for the r-function 



The above simple example already contains all the basic features of at least all the 
Ap cases. Indeed, the reduction ( ^Bj ) is nothing but S'L(2)-reduction of a generic GL{2) 
situation. In other words, the diagonal U{1) G'L(2)-current J = ^{J^^^ + J^"^^) = -^dcj) 
decouples. This is an invariant statement which can be easily generalized to higher p 
cases. 

In the case of SL{p) we have to deal with the p-component hierarchy and instead of 
(|50|) for generic r-function one has 

r^P\x) = (Ar|e^(")G|0) (59) 
N ={iVi,...,iVp}, x = {x^'\...,x^P^} 

and now we have p sets of fermions {tp^'^^*, ^p^^^} i = 1, ■■■,p- The Hamiltonian is given by 
Cartan currents of GL{p) 



H{t) =j:T.^PjiP (60) 

1=1 k>0 



and the element of the Grassmannian in the particular case of CMM is given by an 
exponents of the other currents 



Note that the idea to preserve both indices in (|5^) leads immediately to additional insertions either 
of ■0'^^^* or t/i^^' to the right vacuum |0), so that it is no longer annihilated at least by the r_i Virasoro 
generator, or in other words this ruins the string equation. Thus only the particular reduction (^5|) seems 
to be consistent with string equation. This choice of indices just corresponds to that considered originally 
in M. 
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i.e. 

G = nexp(Q(*^))exp(g(*^))exp(g(*^>) (62) 

gfe) = j gfe') = ^ jfoO^ gfo> = ^ jfe>, z ^ j 

Since (|6l|) are the SL{p)i Kac-Moody currents, (|6^ play the role of screening operators in 
the theory under consideration. It deserves mentioning that they are exactly the SL{p) 
(not GL{p)) -screenings and thus the r-function ( |5UD does not depend on {J2^=iXk^}, i-e. 
we obtain the constraint (|36|) . 

In the case of SL[3) this looks as follows. The screenings are 

Qio^) = j jiOL) ^ (63) 

where {o:} is the set of the six roots of 5*1/(3). In terms of fermions or bosons the screening 
currents look like 







02) 






= tlj^'^>ij^^> = exp(02 + 


03) 




J((X,) 


= ^(1)^(3)* ^ exp(03 - q 


H) 






= ^W^m = exp(-0i - 


02) 


(64) 


J(-(X,) 


= ^/;(2)^(3) = exp(-02 - 


03) 






= ^(3)^(1)* ^ exp(0i - q 


^3) 


(65) 



The particular r-function is now described in terms of the correlator 

tI^\x) = (ATle^^^^^GIO) (66) 

with 

G-Hexpg^"^^ (67) 

The condition of Cartan neutrality is preserved by compensation of charges between the 
operator (|^) and left vacuum (7V| in (|66|) . It is obvious that in such case due to the 
condition of Cartan neutrality of the correlator (like in Wess-Zumino models) the mode 
J = d(j) = Z^iLi decouples from the correlator, and 

r'^^\x) = (Ar|e^(")G|0) = 
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(0|expK:j,$:4M|0)^(iV|e^(*'*"^G|0) 



,fc>0 i=l 



(68) 

= r'N%i), (69) 

where the first correlator in the second row is trivially equal to unity. For the specific 
choice of the operator G in (^) 

G = Gi,2 = exp (^J J(«^)) exp (^J (70) 



we reproduce the formula (^31) . 

Thus, the general statement is that for Ap-series we obtain an "AKNS-type" SL{p)- 
reduction of p-component KP hierarchy. Generalization to other groups is easy. Indeed, 
starting from the corresponding VT-algebra we automatically end up with the element 
of the Grassmannian with proper group properties. Thus, we have a tool to construct 
hierarchies of the matrix model type (i.e. satisfying a string equation) with a given 
symmetry group. All properties of these hierarchies are automatically dictated by the 
symmetry. Say, the defining property of (generalized) AKNS hierarchy is the independence 
of the corresponding p-component KP r-function of the sum of times WM. On the other 



hand, it really corresponds to SL{p) reduction from theoretical group point of view ||25|| . 

It is even more interesting to remark again that the general AKNS system is defined 
for arbitrary set of zero times: 

rN,M ~ (Ar| . . . |M) 

with the only restriction J2 = J2 due to the charge conservation law. But the 
condition that T_i annihilates the right vacuum requires M to be zero (it is interesting to 
note that the original paper [^, devoted to multi-component KP hierarchy, was dealing 
with exactly this type of restriction). In the SL{2) case this gives rise to a one-parameter 
r-function Tjv,-Ar) what allows one to rewrite the system as a Toda (chain)-type hierarchy, 
since the Toda-type system should depend only on a single zero-time (in this case the 
determinant of the type (|i3| ) can be re-ordered to have a block form with the given 
symmetry properties between different blocks). In particular, the first equation from the 
2-component hierarchy which is of the form (|48|) ||29[| , in this case transforms into eq. (^0]) 
which is just a Toda chain equation. 
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As to the Toda-like representation of CMM, in the simplest S'L(2)-case the result is 



indeed equivalent to the Toda chain hierarchy & IT8|, IT9|. In the fermionic language this 



connection is established by the following substitution in the element of the Grassmannian 

i;W(z)^i;(z), (71) 

Z 

and the same for ?/'*'s. This is a reflection of the fact that Toda system is described 
by the two marked points (say, and oo) and corresponds to two glued discs, so it can 
be also described by two different fermions. This might lead to a general phenomenon, 
when any multi-component solution to CMM is actually related to (some reduction) of a 
multi-component Toda lattice. 

4 Double-scaling limit 

In this section we consider the central issue of the connection between matrix models and 
2d gravity theory - continuum limit. Even in the most investigated case of Hermitian one- 
matrix model this is a rather sophisticated procedure, especially if one wants to reproduce 
the whole continuum integrable structure of 0. Moreover, in the case of conventional 
multi-matrix model this procedure is still unknown. 

Below we demonstrate that the advantages of CMM make it possible to define honestly 
the double-scaling limit in these theories along the line of [0 (see also |16[)- This means 



that the set of times of discrete model undergo a Kazakov-like transformation to the con- 
tinuum times. Discrete VT-constraints are transformed into constraints (|I]) of continuous 
model. 

4.1 Results of [7] for the one matrix model 

To begin with let us briefly remind the main points of |]^. 

It has been suggested in that the square root of the partition function of the 
continuum limit of one-matrix model is subjected to the Virasoro constraints 



^cont^/^ = 0, n > -1, (72) 
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where 



^/ Ol2(k+n)+l o<fe<n-l '^-'2A:+lC-'2(n-fc-l)+l 



k=0 



+ 16 + (16G) ^^^^ 
are modes of the stress tensor 

T{z) = --m^z): - ^ = V (74) 

It was shown in that these equations which reflect the l^'^^-'-invariance of the par- 
tition function of the continuum model can be deduced from analogous constraints in 
Hermitian one-matrix model by taking the double-scaling continuum limit. The proce- 
dure (generalized below to CMM) is as follows. 

The partition function of Hermitian one-matrix model can be written in the form 

Z{tk} = IVM exp Tr t^M'' (75) 



and satisfies pO| , PH] the discrete Virasoro constraints (^) 



Ll^Z = 0, n > 



f,^Q Uik+n 0<k<n ^'^k^'^n-k 

In order to obtain the above-mentioned relation between VT-invariance of the discrete 
and continuum models one has to consider a reduction of model (|75]) to the pure even 
potential t2fc+i = 0. 

Let us denote by the rj^'^ the partition function of the reduced matrix model 



red 



A;=0 

and consider the following change of the time variables 



{t2k} = J VMexp Tr J2 hkM^^ (77) 



(-)""'"r (n+ I) fl-"-^ 

9m=Yl — / ^ 1\ ^2n+l, (78) 

n>m [n - my.T [m + 
where Qm = ^t2m, and this expression can be used also for the zero discrete time go = N 
that plays the role of the dimension of matrices in the one-matrix model. Derivatives with 
respect to t2k transform as 

dt2k „t^o(A;-n-l)!r(n + |)9f2„+i' 
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where the auxihary continuum times T2„+i are connected with "true" Kazakov continuum 
times T2n+i via 



T2k+1 — T2k+1 + Q ,^ , 1 /o ^2(fc-l)+l; (80) 



k 

'k + 1/2' 

and coincide with T2„+i in the double-scahng hmit when a ^ 0. 

Let us rescale the partition function of the reduced one-matrix model by exponent of 
quadratic form of the auxihary times T2n+i 

f = exp I -i ^mnf2m+lf2n+l ] t]^'^ (81) 

\ ^ m,n>0 / 

with 

^ ^ r(n+|)r(m+|) (_)n+n^^-n-n^-l 

""^ f n!m!(n + m + l)(n + m + 2) ' 



,2 

Then a direct though tedious calculation [0 demonstrates that the relation 



is valid, where 



r ~ n+l rred^red 

^ = a- 5: C„%,(-1)'^-^^-^^^, (83) 



fc=0 0'r2(fc+n) o<fc<n 0'^2fcCt2{n-fc) 



and 



^-1= E(^+^)^2.-M-^^ + 



2; ^'="^9T2(._i)+i 16G' 



fe>l ^ (>'J2(A:-1) 

^0= E (^+^)^2fe+l-^ 
fe>0 ^ CJ2fc+l 



fc>0 ^ 2^ f^^2{A:+n)+l 

+ ^ i-):, „>i. (85) 

0<fc<n-l dT2k+i dT2 



Here Cf, = ^| are binomial coefficients. 



These Virasoro generators differ from the Virasoro generators (|73|) ^ by terms 

red 
2p 



which are singular in the limit a — > 0. At the same time U2^t'^^'^ at the r.h.s. of (|83D do 



not need to vanish, since 



= L2pT 



i2k + l—0 



T red^rcd , \ ^ 
^2p ^ +2^ 



dt2i+ldt2{n~i-l)+l 
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(86) 



It was shown in that these two origins of difference between (|73|) and (^) actually 
cancel each other, provided eg. (pS]) is rewritten in terms of the square root -\/f rather 
than f itself: 

rcont n+l T ^ T 

p=0 



(l + 0(a)). (87) 

*2fc+l— 



The proof of this cancelation, as given in |^, is not too much simple and makes use of 
integrable equations for r. 

In our consideration of CMM below we will use a more economical way to define 
the change of the time- variables t — > T (also proposed in 0), implied by the scalar 
field formalism. The Kazakov change of the time variables ( [75|J75| ) can be deduced from 
the following prescription. Let us consider the free scalar field with periodic boundary 
conditions ((pHD for p = 2) 

d^iu) = Y: 9ku''~' + E (88) 

fc>0 k>l '^''2fc 

and analogous scalar field with antiperiodic boundary conditions: 



E((* + i)j'.. 

fc>0 \ ^ 



d 



3 



^d^z)f = a-^d^iuy^"^, u^ = l + az (90) 



Then the equation 



generates the correct transformation rules (0), ( ffOD and gives rise to the expression for 
Anm (p2D. Taking the square of the both sides of the identity (pO|), 

lnz)f = :^T(.)r-, (91) 

one can obtain after simple calculations that the same relation (p3|) is valid. 



4.2 On the proper basis for CMM 

In this subsection we would like to discuss briefly the manifest expressions for constraint 
algebras of the sect. 2 in terms of time- variables. Indeed, for convenience of taking the 
continuum limit, the time variables should be redefined {i.e. the integrable flows of the 
previous section are not suitable in the continuum limit). In other words, this is the 
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question what is the proper reduction, or what combinations of the "integrable" times 
should be ehminated. 

To begin with, we consider the simplest non-trivial case of p = 3. Then introducing 
the scalar fields 



)=j:kt','z'-'+j:^^-'-\ (92) 



k k Ot^ 



k k Otf. 

with t^f}^ = {itk + tk)/2y/2, t^^^ = [itk — tk)/2^/2, one obtains the expressions: 

iy(2)(^) = ]^d(P^^\z)d(t)^''\z), (94) 

W^^'n^) = ^E(^</>^^^W)'- (95) 

instead of (H) and (H). 

This choice of basis in the Cartan plane is adequate to the continuum limit of the 
system under consideration, as the latter one is described by completely analogous expres- 
sions 0. Now let us describe this basis in more invariant terms and find the generalization 
to arbitrary p. 

Comparing (P^) with (EH), we can conclude that S^'-*-' = f3^d(f) corresponds to the basis 



/3i,2 = ^(V3M2±^«2). (96) 

This basis has the properties 

/3i-/32 = l, /3i-/3i = 0, /32-/32 = 0. (97) 

Now it is rather evident how this basis should look in the case of general p. Due to |0] 
we can guess what is the choice of the proper scalar fields: 

k k Ot^. 

This choice certainly corresponds to the basis with defining property (it can be observed 



immediately from the relations (|98|) and (|18])): 

/3, ■ /3^. = (99) 
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the proper choice of the Hamiltonians in (|T7D being 

H = J2t^f3rJk, (100) 

i,k 



what determines new times adequate to the continuum hmit. 

Let us construct the basis ( p9D in a manifest way. To begin with, we define a set of 
vectors {/Xj} with the property: 

Hi-Hj = 6ij--, ^/ii = 0. (101) 

P i 

The H^(")-algebra can be written in this basis as follows B: 

n 

^(n)^(_)n+i ^ n(%„-50), n=l,2,...,p. (102) 

l<ji<...<jfc<p m=l 

Now the basis ( pQ]) can be constructed from ( |101j ) by diagonalization of the following 



cyclic permutation g, |33 



Mi^/^i+i> Mp^Mi i = l,...,p-l. (103) 

This transformation has as its eigenvectors, their manifest expressions being of the 
form: 

1 ^ 2TTi 

/3, = — ^ exp{—jk}fx^, A; = 1, 2, . . . , p - 1. (104) 

It is trivial to check that the properties (|99D are indeed satisfied. One can immediately 
rewrite the corresponding ly-generators in the basis of /3j's. After all, one obtain the 
expressions similar to the continuum W^-generators 0, but with the scalar fields 
defined as in ( pSD and without the "anomaly" corrections appearing in the continuum case 
due to the twisted boundary conditions. These corrections can be correctly reproduced by 
taking the p-th root of the partition function as well as simultaneously doing the reduction 
(see subsects.4.3-4.5). 

Thus, the proposed procedure allows one to take the continuum limit immediately 
transforming the scalar fields as elementary building blocks. 



4.3 The two- matrix model example 

In this section we consider the simplest example of the p = 
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3 CMM. 



From {t}- to {T} -variables. To describe this change of variables we shall use the 
scalar-field formalism. 

Consider the set of scalar fields (p^) , (p3D and perform the reduction 

tg+i = 42+2 = 0, ^ = 1, 2, A; = 0, 1, . . . . (105) 

Then 



fc>0 fc>l C'^sfc 



d^'^'in) = E 9^^'-' + E i^"''"' (106) 

A:>0 k>l Ut^^ 



and 5'^*'' = Hgj, = 1,2,..., gfg*^ = A^*^*-* were introduced in the subsect.4.2. A^(^'^) are 
"the dimensions of the matrices" used in our two-matrix model. Then we put 

^9$^ {z)f = au'-^^dip^'^ {uy"^, u^ = l + az, i = 1, 2, (107) 

where ^^^Hz) are the scalar fields of the continuum model 01 



fc>0 V ^ CJ3A:+2 



and 

k 

'k + i/?, 



T^k+i = Tsk+i + a^— p— 777T'3(fc_i)+i, i = 1, 2. (109) 



The relation ( |107|) gives rise to the following Kazakov-like change of the time variables 



9^^ = E r3n+. ^ = 1, 2 (110) 



n>m (n-m)!r(m+|) 



and equations ( p.lO| ) can be also continued to the zero discrete times gg K It follows from 
(|107|) that the derivatives with respect to t^^l, i = 1,2 transform as 



= !^ U Jl i = 1 2 (111) 

dt'il t'o{k-n-ly.T[n + l + l)^nn+^' 

To cancel the terms like T3k+2Z~^~^ and T^k+iz^''^^ with /c > at the left hand sides of 
( |107| ) we have to rescale the partition function of the continuum model f by means of an 
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exponent of some quadratic form 



f = exp|- Yl ^mnT3^+2T3„+l I r'^^^ (112) 
m,n>0 



where A„m has the form 



\n+m ^—n—m—l 



rfn + f)r(m + |) , , ^ 

r(^)r(i) n\m\{n + m + l){n + m + 2)' ^ ' 



Constraint algebras of the two-matrix model. It was suggested in that the continuum 
"two-matrix" model possesses the W*-'^^ (including Virasoro) symmetry, where Virasoro 
generators and generators of the W^^-* algebra are constructed from the scalar fields 
(|108|) in the following way 

nz) = l:9$(^)(.)9$(^)(.): - ^ = E (114) 

W^'\z) = ^ (: (9$«(.))' : + : (5$(^)(.))^ ^ = ^ ^^^^^ 

One can easily derive a relation between the generators and the corresponding 
generators L^^ny associated with reduction ( |105| ): 

-| n+l rred—red 

= a-"ECr'(-r'''"'^^' ^>-l' (116) 

where Virasoro generators are defined by the same formula (p^), only with substi- 
tuted by "reduced" fields ip: 

Tin) = ^:av.«(n)V')(«): = E«"'"~'^3;^ (117) 



and 



EI(* + 5)j-»«5£;H-(*4)^»«5Tb) -^5^'^- 



fc>l 



fc>o \^ "^-^^ fW^3fe+i ^ 3/ dTsk+2/ 



+ E — — ^ + ^^1- (118) 

m+fc=n-l dT^m+l 9^" 
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Eq. (|116|) is a direct consequence of the relation 

=l + az (119) 



_rcd 



(compare with (|107| )). As in the 1-matrix case, the generators ( |118| ) differ from the 



continuum generators of by singular c-number terms. Instead, again do not 
exactly annihilate r'''^'^. We assume that again these two effects cancel each other, provided 
eq. ( |116| ) for f is rewritten in terms of the cubic root v^- In orther words, doing accurately 
the reduction procedure in the Virasoro constraints of the discrete two-matrix model one 
can rewrite (|116|) in the form 



1 "^^^ T T 

^/:r*v^ = a"" E C;+^(-)"+i-^' ^ (1 + 0{a)) 
y/T p=Q T 



, n > -1, (120) 

i3fc+i=0; i=l,2 



Thus, we conclude that the Virasoro constraints of the continuum matrix model are indeed 
implied by the corresponding Virasoro constraints of the discrete conformal two-matrix 
model. 

It follows from ( |115|) that generators of the W^^^ symmetry for the continuous model 
can be written in the form 

nX3) = 3 E (k + I) (m + 1) r3,+iT3„+i ^ 



k,m>0 ^ 3J \ 3/ 9T3(fc+m+n)+l 

+ 3 E (m+'^)T3k+2T3m+2—~ 

V- / 4\ d d 
+ 3 2^ (m + p-n + -]T^m+p-n+i)+i^ 

m,p>0 ^ C'-' 3m+2 C-f 3p+2 

V- f 2\ d d 

+ 3 2^ [m+p-n + -] T3(^m+p-n)+2 



m,p>0 ^ CJSm+lCJsp+i 

+ E S i ^ +^ i ^ , (121) 

n—k—m)— (n— fc— m)— 2 

where n > —2 and the terms with negative values of the indices should be omitted. The 
generators VV_2 and Wi/ (but not the Wq ) have additional terms, cubic in times: 

W^^l = + . . . and W^^l = + ^TlT^ + ■■■ . (122) 
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Similar to the Virasoro case one can show that the relation between the generators of 
the W^^^ symmetry for the discrete model and the generators VV^'^^ is 



1-r 



P+2 W^^^ red red 



n=0 



r 



red 



p>2 



(123) 



and follows from the identity 



T 



^3 ^ T-red 



UjT 



red ^ 



\ -j-red 



red ^ 



(124) 



where the generators PV^^-* '^^'^ of the discrete model are defined after the reduction ( |105| ) 
by the relation 

1 



-3n-3u/(3) red 
3n 



(125) 



After reduction of the discrete model eq.( |123| ) can be rewritten in the form 



1 P+2 



P>2, 



(126) 



t3fe+,=0; i=l,2 



where vy^^-*™"* are the generators of the W^^^-symmetry from the paper (see ( |121D and 
f|T22D). 



4.4 The general case 

It is easy to generalize the example of p = 3 to general multi-matrix models using the 
formalism of scalar fields with Zp-twisted boundary conditions. Let us introduce p—1 sets 
of the discrete times t^*"*, i = 1, 2, . . . , p — 1 and A; = 0, 1, ... for the discrete {p — l)-matrix 
model and consider the reduction 



''pk+j 



0, ^,j = l,2,...,p-l, k = 0,l, 



(127) 



Choose the discrete and continuum scalar fields in the form 



d 



k>0 



k>l Otp,^ 



-pk-1 (i) 



iJk — '^''pk 



d 



a$«(.) = EU^ + - 



Z V 



k>0 



(128) 
(129) 
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aM'-P+^^9(^» {uy"^ = {z)f, uP = l + az, i = 1,2, . . . ,p - 1, (131) 



Tpk+i = Tpk+i + aj—^—rj^Tp(k-i)+i, i = 1,2, . . . ,p — 1. (130) 
Then the equations 

l^V^)(«)r-'^ = l 

generate the Kazakov-hke change of the time variables 

(-)"-"T fn + 1 + ^) a~""p 

9^= E : ^ = l,2,...,p-l (132) 

n>m {n - my.T [m + -j 

-^ = y ^ H A— , i = l,2,...,p- 1. 133 

dt^;l to{k-n-l)\T{n+l + l)dTp^+- 

Using the eq. (|131|) and considerations similar to those used in the previous subsections 
one can show that there is relation between tilded continuum generators VV^*-*, i = 2, . . . ,p 
of the W-symmetry and reduced generators of the discrete VT-symmetry W^^'^^'^ of the 
form 

1 _ n+j-l Tl/(i)red^red 

= a-" Cr-^(-)"+-^- , n>-z + l, (134) 

where rescaled r-function is defined 

f = exp I -if: E Al^nTpm+^Tpn+p^i] T^'^^itpk} (135) 
y i=l m,n>0 J 

and matrices A^*^ are determined by 

^"'"^ r(i)r(2^) n\m\{n + m + l){n + m + 2)' ^ = 1' 2, . . . ,p- 1. (136) 

This relation again corresponds to the identity 



(137) 



where = 1 + az. 



Performing the proper reduction procedure (|127|) , which eliminates all but the time- 
variables of the form tpl (i.e. leaves the 1/p fraction of the entire quantity of variables) 
we can obtain the relation 



n+i-l W^^^T 
-n+j-1/ \n+i~l-s " ps ' 



ii^O only 



(l + 0(a)), n>-t+l, 

(138) 



29 



where W^*-' is the W-generators of the paper . Thus, we proved the 14^-invariance of the 
partition function of the continuum p — 1-matrix model and found the exphcit relation 
between its partition function and corresponding partition function of the discrete (p — 1)- 
matrix model. 



4.5 On the reduction of the partition function 

To conclude this section we would like to discuss the problem of reduction 



and 



(|127|) of the partition function in detail, with accuracy up to (non-leading) c-number 
contributions and only after the continuum limit is taken. More precisely, we reformulate 
the condition of a proper reduction in the continuum limit in order to reduce it to more 
explicit formulas which can be immediately checked. As a by-product of our consideration 
we obtain some restrictions on the integration contour in the partition function (p2D. 

To get some insight, let us consider the simplest case of the Virasoro constrained 
Hermitian one-matrix model M. Before the reduction the Virasoro operators read as in 



([76|). Then their action on logr can be rewritten as 



(9 logr 



E 



(9^ log r 



After the reduction, we obtain 



+ E 



91ogr'-'=<i ^ 92 log r--^^ 



L k 



+E 



dt2k+2n m 9t2mdt2n —2m 



+ E 



d log T d log r 
~~df 'df 

Q2 logr'-^d 



(139) 



dt2m+ldt 



2n-2m~l 



+ 



dt 
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under the condition 
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The last formula is a direct consequence of the "Schwinger-Dyson" equation induced 
by the transformation of the reflection M — > —M in ([751). Indeed, due to the invariance 
of the integration measure under this transformation one can conclude that the partition 
function ([75|) depends only on a quadratic form of odd times. 

Thus, the second derivatives of logr over odd times do not vanish, and are conjectured 
to satisfy the relation 

g2 ^red 



(9^ log- T^^'^ 
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dt2m+ldt2n-2m-l 
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where the sign ~ imphes that this relation should be correct only after taking the contin- 
uum limit. In this case one obtains the final result (cf. (p^)) 



■ d log Vt^ ^ 92 log Vt^ 



0. (143) 



dt2k+2n m 9t2mdt2n-2m _ 

Thus, it remains to check the correctness of the relation ( p.42| ). To do this, one should 
use the manifest equations of integrable (Toda chain) hierarchy, and after direct but 
tedious calulations [0 one obtains the result different from the relation ( p.42| ) by c-number 
terms which are singular in the limit a — and just cancell corresponding items in (|85D 
(this is certainly correct only after taking the continuum limit). 

All this (rather rough) consideration can be easily generalized to the p-matrix model 
case. In this case one should try to use all W^^-'-constraints with 2 < i < p. Thus, 
the second derivatives should be replaced by higher order derivatives, and one obtain a 
series of equations like ( |139| ). It is the matter of trivial calculation to check that these 
equations really give rise to the proper constraints satisfied by ^^r (cf. (|120|) , (0) and 
(|138|) ) provided by the two sets of the relations like ( |141|) and ( |142| ) . 



Namely, the analog of the relation ( |141| ) in the p-matrix model case is the cancellation 
of all derivatives with incorrect gradation, i.e. with the gradation non-equal to zero by 
modulo p. The other relation ( |142|) should be replaced now by the conditions of the 



equality (in the continuum limit) of all possible terms with the same correct gradation. 
In the simplest case of p = 3 these are 



log r'"'^'^ ff^ log 
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m.k 9t3rn+2dtsn-3{m+k)+ldt3k ^^j^ dt3rndt^n-Z{m+k)dt2,k 

Again, this second condition is correct modulo some singular in the limit of a — > 
terms, which appear only in the case of even p. Unfortunately, we do not know the way to 
prove this statement without using the integrable equations, what is very hard to proceed 
in the case of higher p. 

On the other hand, the cancellation of derivatives with incorrect gradation can be 
trivially derived from the "Schwinger-Dyson" equations given rise by the transformations 
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M exp |2^|M (0 < k < p) of the integration variable in the corresponding matrix 
integral, the integration measure being assumed to be invariant. In its turn, it implies that 
the integration contour, instead of real line, should be chosen as a set of rays beginning 
in the origin of the co-ordinate system with the angles between them being integer times 
This rather fancy choice of the integration contour is certainly necessary to preserve 
Zp-invariance of p-matrix model system. 

5 Conclusion 

In this paper we proposed a new point of view on the discrete matrix formulation of the 
unified theory of 2d gravity coupled to minimal series of c < 1 matter. As alternative 
to conventional discrete multi-matrix models which are extremely hard to investigate by 
imposing some differential equations (Ward identities etc.) on matrix path integral, we 
introduced another class of models which satisfy "conformal" constraints by definition and 
actually possess much richer integrable structure, essentially as rich as the one arising in 
continuum limit. This is also up to now the only case when continuum double scaling 
interpolating limit (of multi-matrix models) can be performed honestly, reproducing the 
results of 0. 

If matrix models are considered as certain solutions to integrable theories, the confor- 
mal multi-matrix models satisfy the equations of mult i- component KP hierarchy. More 
precisely, they correspond to particular reduction of multi-component KP hierarchy, the 
reduction of AKNS type, which is further constrained to be consistent with discrete W- 
constraints and/or discrete string equation. We found a determinant representation for 
a certain subclass of such solutions, which can be also considered as a generalization of 
those corresponding to orthogonal polynomials, leading to forced hierarchies. All this 
must shed light on the origin of the integrability in 2d gravity and the work only started 
in this direction. 

So far there was no real progress in taking the continuum limit of conventional multi- 
matrix models, and most of the facts were rather introduced axiomatically in this case. 
In contrast to this, CMM do have a nice continuum limit, which can be described in the 
same terms as the continuum limit of 1-matrix case. This shows that the fact of discrete 



32 



ly-invariance or, put differently, the proper choice of a particular representative within 
the universality class, is crucial to expose the origin of string equation of 2d gravity theory 
in a manifest form. 

Certainly there must be a deep connection between CMM and unified description of 
the continuum theory via Generalized Kontsevich Model (GKM) |^ |^, |T9|, ^ |3^. We 
know 1 19 1 that there exists also a reformulation of a discrete Hermitian 1-matrix model in 



terms of GKM, the integrability being established using formulas very close to those of the 
sect. 3. It is possible also to take the continuum limit of 1-matrix case in "internal GKM 
terms", what is much easier than by technique developed in |0] and sect. 4 above. All 
these facts implies that there should be a sort of reformulation of the CMM in GKM-like 
terms. 

There is another still unresolved question on interpretation of conformal fields ap- 
pearing in the above formalism and their relations to the conformal fields of Polyakov's 
formulation of 2d gravity. Of course, there shouldn't be any direct correspondence, but 
there can be a kind of duality between "world-sheet" and "spectral" Riemann surfaces. 
The key role of 14^-constraints found in the paper and the presence of multi-component 
scalar fields might also have a meaning in the framework of 1^-gravity theory |3^, |3^ . 

We hope to return to all these problems elsewhere. 
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